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We consider the simplest versions of the Nambu-Jona-Lasinio (NJL) model and the Linear Sigma 
Model (LSM), in the Mean Field Approximation (MFA), in order to analyze hot and dense two 
flavor quark matter subject to strong magnetic fields. We pay especial attention to the case of a 
finite chemical potential, which has not yet been fully explored. Our results, for the NJL model, 
are in qualitative agreement with other recent applications showing that, for stronger fields, the 
first order segment of the transition line increases with the magnetic strength while the coexistence 
chemical potential value, at low temperatures, decreases. In the present work, one of the most 
- - - important results is related to the analysis of how these features affect the phase coexistence region 

(N- in the T — pB plane. We find that the coexistence boundary oscillates around the B = value 

, for magnetic fields of the order eB < 9.5 which can be understood by investigating the filling 

>— ^ ■ of Landau levels at vanishing temperature. So far, most investigations have been concerned with 

the effects of the magnetic field over the T — fi plane only while other thermodynamical quantities 
{3JQ[ such as the adiabats, the quark number susceptibility, the interaction measure and the latent heat 

^ ■ have been neglected. Here, we take a step towards filling this gap by investigating the infiuence of 

a magnetic field over these quantities. Finally, we argue that a naive application of the MFA does 
not seem to be appropriate to treat the LSM in the presence of magnetic fields. 
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O : I. INTRODUCTION 

ay 



The determination of the QCD phase diagram, even in the absence of magnetic fields, is still a matter of great 
theoretical and experimental activities. In this case, powerful lattice simulations have established that at vanishing 
C^l I baryon chemical potential there is no true phase transition from hadronic matter to a quark gluon-plasma but rather 
a very rapid raise in the energy density signaling a crossover characterized by a pseudocritical temperature, Tpc ~ 



' 160 MeV The situation is less clear for the finite chemical potential region since, so far, there is no reliable 
. information available from lattice QCD evaluations. Nevertheless, most finite /i lattice extrapolations for the /i = 

■ Columbia plot indicate that the critical first order surface (in the mu^d — n^s — M space) will hit the physical current 
mass values at some finite ^ thereby characterizing a critical end point (CP) Q which constitutes the most plausible 
theoretical scenario. On the other hand, other lattice evaluations predict that, as increases, this critical first 

[ order surface bends in a way which in principle would exclude the appearance of a CP and the whole T — ^ plane 
would then be dominated by a crossover representing an "exotic" scenario [2]. However, even this situation may be 
reversed, in favor of a critical end point, if a different physics strongly influences the finite density region causing a 
back bending of the critical surface S-'g] . 

■ At the same time, the possibility that strong magnetic fields may be produced in non central heavy ion collisions 
Hi as well as being present in magnetars iSj , and in the early universe leads to the question of how these fields 

influence the QCD phase diagram. So far, most estimates have been carried out at vanishing chemical potential with 
the aid of effective theories such as the Linear Sigma Model (LSM) [l^] and the Nambu-Jona-Lasinio model (NJL) 
[nilll within the Mean Field Approximation (MFA). The general outcome is that the pseudocritical temperature, 
at which the crossover takes place, increases with increasing values of the magnetic field. Some of these applications 
have also considered these effective models with the inclusion of the Polyakov loop to take confinement into account. 
The results show that the Tpc related to deconfinement also increases with B but an interesting splitting between 
this temperature and the one related to chiral symmetry has been observed At /j, = 0, the first lattice attempt 
to solve the problem considered two quark flavors and high values of pion masses (m^r = 200 — 400 MeV) confirming 
that Tpc should increase with B [l^. However, an improved lattice simulation [T4| which considered 2+1 quark 
flavors at physical pion mass values (m^ = 140 MeV), together with an extrapolation to the continuum, predicted 
that Tpc should decrease with B. A decrease of Tpc, for the deconfinement transition, has also been observed with 
the MIT bag model fld\ and, more recently, with the large- TVc limit of QCD [16j| (see Ref. 17 1 for a summary of these 
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results) . In Ref. [la] it is suggested that this behavior is due to the entanglement with the Polyakov loop that is 
affected by magnetic screening as discussed in Ref. [Toj . where it was shown, within a simplified framework, that the 
strong electromagnetic fields can play a catalysing role for a deconfinement transition. However, as referred in Ref. 
(Ts"], if only the chiral sector is considered, an enhancement of the chiral condensate, due to finite B, would not be 
contradictory with the lattice simulations of Ref. [141 . Nevertheless, although the quantitative discrepancy remains 
to be fully understood, there is a concensus regarding the crossover character of the chiral transition at vanishing 
chemical potential. On the other hand, the effects of strong magnetic fields at the finite fi regime has not been fully 
explored to date even within effective theories. Experimentally, the intermediate temperature regime is important 
regarding current nuclear collisions at the low-energy end of RHIC and in the future with FAIR at GSI and NICA at 
JINR whose aim is to probe the expected critical point region. At the low temperature end one expects to find a rather 
interesting physics related to the expected first order phase transition which may have astrophysical consequences 
(e.g., regarding the possibility of quark star formation |20l. [2]|). 

In the chiral limit, the two flavor version of the NJL model subject to a magnetic field was considered in Ref. [l^ 
with the MFA while the three flavor version, at the physical point, was recently analyzed in Ref. [l^ within the same 
approximation. In the latter reference, which concerns the more realistic case, it was observed that Tpc, at /i = 0, 
always increases with B, as in most model applications. In that work, one of the main novelties regards the critical 
end point location, which moves towards higher temperatures and smaller chemical potential values as B increases. 
Another interesting result concerns the size of the first order segment of the transition line, which becomes longer, as 
stronger magnetic fields are considered. Also, at low temperatures, it has been observed that the chemical potential 
value at which the first order transition occurs decreases with B. The latter result has been previously observed with 
the two fiavor NJL, in the chiral limit [l^, as well as with a holographic one-flavor model [Ij] . The interested reader 
may find a model-independent physical explanation for this result, which the authors have termed Inverse Magnetic 
Catalysis (IMC), in Ref. [l^] while a recent review with new analytical results for the NJL can be found in Ref. f25| . 
Also, very recently, the two flavor NJL model has been considered to investigate the dynamics of neutron mesons in 
a hot and magnetized medium [2^. Another important application regards the generalized three flavor NJL model 
including eight-quark interactions in which the phenomenon of secondary magnetic catalysis sets in . 

The LSM with two flavors has also been employed to determine the T — /i phase diag ram with the MFA [2^ 
and, more recently, with the more powerful Functional Renormalization Group (FRG) [29| . This method has been 
previously used in an application to the Polyakov quark model (PQM) demonstrating that, at fj, = 0, the increase 
of Tpc with B persists even when mesonic fluctuations are considered [3(i]. The MFA results of Ref. [2^ show that, 
when all fermionic contributions are considered, the whole T — fi plane is dominated by the crossover, irrespective of 
the magnetic field value while the chiral symmetry broken region expands with B. It is interesting to note here that, 
in this case, the LSM mimics the "exotic" theoretical QCD phase diagram scenario mentioned at the beginning of 
this section. However, as we show here, this is just an artifact of the parametrization adopted in Ref. p8| where the 
authors chose a rather high value for the sigma meson mass (m^ — 800 MeV). As a matter of fact, the occurrence 
of a critical point and first order transition depends crucially on the mass of the sigma meson that one uses in 
the computation [sT'l. Apart from the parametrization issue one must be very careful when considering the LSM 
since different approximations may lead to very different results even at the qualitative level (e.g., yielding different 
types of phase transitions 0, H^, H^l). This is mainly due to the fact that within this effective theory one has 
scalar and pseudoscalar mesons as well as quarks as degrees of freedom. Also, the pure vacuum contributions may 
be neglected without spoiling the breaking of chiral symmetry, which happens at the classical level. This allows 
for different approaches in which one or more contributions are neglected yielding different results. When magnetic 
fields are present, the phase diagram obtained with the FRG applied to the LSM j29|] agrees, qualitatively, with the 
ones obtained with the NJL in the MFA [2^, [2^. However, sigma masses of the order nia = 400 — 450 MeV have 
been considered so that one may wonder if the qualitative agreement between the LSM-FRG and the NJL-MFA 
predictions for the chiral transition phase diagram are due to the use of smaller rria values or to the use of a more 
powerful approximation scheme within the LSM. Here, in order to address this question, we consider the LSM with 
the MFA using mo- = 600 MeV and also mo- — 450 MeV showing that the higher value reproduces the "exotic" scenario 
also found in Ref. On the other hand, using mo = 450 MeV, we observe the appearance of a first order line 

starting at T = and terminating at the critical end point. However, when magnetic fields are turned on, we do 
not observe the same qualitative features observed with the NJL-MFA and the LSM-FRG at intermediate to low 
values of the temperature. For instance, at T = 0, the coexistence chemical potential values does not decrease with 
B as predicted by the Inverse Magnetic Catalysis mechanism [23, which points out to the fact that the naive 



MFA application as performed here, and also in Ref. [28[, may not be adequate to treat the LSM in the presence of 
magnetic fields. 

Finally, note that despite all the progress made so far in analyzing the infiuence of B over the T — ^ phase diagram 
not much effort has been devoted to investigate particular aspects such as the coexistence and spinodal boundaries 
associated with the first order transition. Also, the analysis of quantities such as the quark number susceptibility. 
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adiabats, interaction measure, latent heat, etc may help in the understanding of the phase diagram structure under 
strong magnetic fields. Therefore, the main goal of the present work is to improve over the simple T — fi phase 
diagram by investigating its structure and related physical quantities in more detail. For instance, to the best of 
our knowledge, the influence of B over the coexistence and the critical regions, for example, has not been addressed 
before. 

In this work we consider the LSM and the N JL model with two flavors in the framework of the MFA in an application 
which could be viewed as an extension of Ref |33] to the B ^ case (the thermdoynamics of both models, at _B = 0, 
has been recently discussed with great detail in Ref. [1^). Note that, here, our main focus is on the NJL which 
yields the expected phase diagram scenario, even at the MFA level, allowing us to probe the physically rich first order 
transition region with its associated critical end point. On the other hand, as already mentioned, the MFA treatment 
of the LSM at B ^ does not generate the expected phase diagram in the full T — /i plane by excluding the IMC 
phenomenon for example. Here, its investigation is justified for allowing us to trace this problem as being generated 
by the MFA and not only by a parametrization which considers high ma- values. These two models, extended by 
the Polyakov loop, have recently been considered in Ref. [sHi where quantities like the magnetic susceptibility of the 
quark condensate as well as the quark polarization have been evaluate at vanishing temperature. 

In the next section we investigate the phase structure of the NJL paying special attention to the low temperature 
regime which has been less explored in the literature. In the same section we evaluate the effect of B over the 
coexistence and spinodal regions, isentropic lines and the quark susceptibility with the aid of the T — fi, T — pB and 
P — T planes. In Section 3, we consider the LSM in the MFA framework, as in Ref. 28], but extending the analysis 
so as to consider the effects of B in the sigma meson mass. We also present the T — /i phase diagram for this model, 
at the physical point, for the magnetic field values relevant to RHIC and the LHC. This exercise will allow us to 
compare the MFA results furnished by both models. Our conclusions are presented in Section 4. 



II. THERMODYNAMICS OF THE NJL UNDER A MAGNETIC FIELD 



At vanishing temperature and finite density the NJL model, subject to a magnetic field, has been used to address 
different questions such as the stability of quark droplets [s^l and the EoS of magnetars with and without strangeness 
in Refs. "21^ and [20| respectively. The influence of B and instantons at the two extreme parts of the phase diagram 
{T — 0, fi ^ and T ^ 0, ji = 0) was studied in Ref. [11]. The whole T — fi plane was first analyzed in Ref. 22] where 
only the chiral limit, for the two flavor version, was considered. At the physical point, the same version of the model 
has been recently considered in Ref. f26l] whose results, for the phase diagram, qualitatively agree with the ones found 
in Ref. [23j for three flavors. In this Section we will obtain the EoS and then analyze different physical quantities in 
order to characterize the physical situations which are more sensitive to magnetic effects. 



A. The Model 

The NJL model is described by a Lagrangian density for fermionic flelds given by [3£ 



W =i'{i^~m)^p + G [(Hj)^ ~ (V^TsTi^)'] , (2.1) 

where -0 (a sum over flavors and color degrees of freedom is implicit) represents a flavor iso-doublet (u, d type of 
quarks) iVc-plet quark fields, while r are isospin Pauli matrices. The Lagrangian density (|2.ip is invariant under 
(global) U{2)f X SU{Nc) and, when m = 0, the theory is also invariant under chiral U{2)l x C/(2)/j. 

Due to the quadratic fermionic interaction, the theory is nonrenormalizable in 3+1 dimensions [G has dimensions of 
eV-2), meaning that divergences appearing at successive perturbative orders cannot be all eliminated by a consistent 
redefinition of the original model parameters (fields, masses, and couplings). The renormalizability issue arises during 
the evaluation of momentum integrals associated with loop Feynman graphs in a perturbative expansion and, in 
the process, one usually employs regularization prescriptions (e.g. dimensional regularization, sharp cut-off, etc) to 
formally isolate divergences. However, the procedure introduces arbitrary parameters with dimensions of energy that 
do not appear in the original Lagrangian density. Within the NJL model a sharp cut off (A) is preferred and since 
the model is nonrenormalizable, one has to fix A to a value related to the physical spectrum under investigation. This 
strategy turns the 3+1 NJL model into an effective model, where A is treated as a parameter. The phenomenological 
values of quantities such as the pion mass {m^r)^ the pion decay constant (/,r), and the quark condensate ({ipip}) are 
used to fix G, A, and m. Here, we adopt the values m = 6 MeV, A — 590 MeV and GA^ = 2.435 which have also been 
employed in Ref. [11] (see Refs. [i3,|4l[ for other possibilities). 
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B. The NJL free energy in the presence of a magnetic field 

The free energy, in the MFA, can be written as fohows [13, (see Ref. m| for results beyond MFA) 

where M is the constituent quarks mass. In order to study the effect of a magnetic field in the chiral transition at 
finite temperature and chemical potential a dimensional reduction is induced via the following replacements in Eq. 



Po i{u>^ ~ ifi) , (2.3) 
pI + {2n + 1 - s)\qf\B , with s = ±1 , n = 0, 1, 2... , (2.4) 



d^P .^Mf^gp^^ (23) 



^ (27r)4 27r ^ ^ J ^ 2Tr 



where = {21/ + l)7rT, with i/ = 0,±1,±2... represents the Matsubara frequencies for fermions, n represents the 
Landau levels and is the absolute value of t he q uark electric charge {\qu\ — 2e/3, \qd\ — e/3 with e = l/-\/137 
representing the electron charge) . Following Ref. |20| we can write the free energy as 



T-NJL _ (-M" to) t-NJL I x-NJL I x-NJL (n (^\ 

' ''vac ' *^mag ' mcd ' \^'^) 



where 



Tli^ -2N^N, I -^{P' + M^fl^ . (2.7) 



This divergent integral is regularized by a sharp cut-off, A, yielding 



jM^ln 


[(A + ea)] 




M 



-eAA[A2 + 4]| , (2., 



vac o 2 

where we have defined ea — VA^ + M^. The magnetic and the in-medium terms are respectively given by 



•^.n^g = jlmBf |c^^^°)(-l,x,) - i[4 - ln(x,) -f ^1 , (2.9) 



and 



j-NJL^_^y-y-^^| |5 I °°^/rin[i + e'[^..'=(^)+'']/^]+rin[l-fe-[^-'=(^)-''l/^]| . (2.10) 

Note that in the last equation we have replaced the label by fc in the Landau levels in order to account for the 
degeneracy factor = 2 — (5ofe. In Eq. (|2.7p . N^. = and Nf = 2 are the color and flavor degrees of freedom, 
respectively. Also, in Eq (|2.9p we have used Xf = M"^ / {2\qf\B) and C,^^'^\—l,Xf) = dC,{z, Xf)/dz\z=-i with ({z, Xf) 
representing the Riemann-Hurwitz function (the details of the manipulations leading to the equations above can be 
found in the appendix of Ref. Finally, in Eq. (PHH) we have Ep^kiB) = ^pl +2fc|q/|B + Af2 where M is the 

effective self consistent quark mass 
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M 



m ■ 



N^NfMG 



d 



AVA2 + M2- ^In 



(A + + Af2)2 



M2 



E |lii[r(^/)] - ^ ln(2^) + - 1 (2a;/ - 1) ln(x/)| 

j: 7 n — C 



N,MG 



-OO -^P,i 



(2.11) 




FIG. 1: Phase diagram for the NJL model in the T — fi plane. The CP, represented by solid symbols, occurs at (Tc ~ 
78.5MeV,Atc ~ 328 MeV) for B ^ 0, {Tc ~ 81.7MeV,^c ^ 318 MeV) for eB = 6ml, and {n ^ 115.6 MeV,/ic ^ 279.3 MeV) 
for eB — 15 ml. The continuous lines represent first order phase transitions and the dashed lines represent crossovers. The 
inset shows the crossing of the transition lines. 




FIG. 2: Left panel: Spinodal boundaries for the NJL model in the T — fi plane. Right panel: Phase coexistence boundaries in 
the T — pB plane (ps appears in units of the nuclear matter density, po = 0.17fm~'^). The solid symbols indicate the location 
of the critical point for each value of B. 



C. Numerical Results for the NJL model 



We can now obtain the phase diagram as well as other important physical quantities for the NJL model in the 
presence of B. All the relevant thermodynamical quantities can be readily obtained by recalling that the free energy, 
evaluated at the mass value which satisfies the gap equation, gives the negative of the pressure, J-{M) = —P. Then, 
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the net quark number density is obtained from p = dP/dfi, and the entropy density from s = dP/dT while the 
energy density is £ = —P + Ts + fj,p. Let us start with the phase diagram, in the T — p, as shown by Fig. [T] which 
was obtained for vanishing B as weh as for eB — 6m^ and eB = 15mJ, which cover the estimated values for non 
central collisions at RHIC and the LHC, respectively [S]. At fi — one observes an increase of Tpc which was 
expected to happen within the present model approximation, as discussed in the introduction. However, at around 
H « 285 MeV the Tpc value for eB = 6 is smaller than the B = value, in accordance with the findings of Refs. 
[23I m, [13 . The figure also indicates that B induces a noticeable increase of the first order segment of the transition 
line which terminates at higher values of T and smaller values of fi. At low temperatures, the coexistence values 
of the chemical potential decrease with increasing B showing IMC (although /z starts to increase again for values 
eB > 16m^ ~ [l^, which we do not consider here). At T = this pattern was also observed in Ref. [37| . 
The CP is located at (T^ ~ 78.5 MeV, - 328 MeV) for B = 0, (T^ ~ 81.7 MeV, /z^ ~ 318 MeV) for eB = 6ml, and 
(Tc ~ 115.6 MeV, /^c ^ 279.3 MeV) for eB = 15 m^. 




FIG. 3: Left panel: The NJL model effective quark mass (squares) at the lowest value occurring at the transition, , and 
the highest coexisting baryon density (dots), pg (in units of po), as functions of eB/m% aXT — Q. The lines are shown just in 
order to guide the eye. Right panel: The latent heat, Af, at T = 0, as a function of eB/ml. As expected, this quantity and 
Pg behave in a similar way. 



100 




T [MeV] 



FIG. 4: The NJL model phase diagram in the P — T plane indicating the regions of broken, CSB, and (partially) restored 
chiral symmetry, GSR, which correspond to the "gas" and "liquid" phase respectively. All lines represent first order phase 
transitions which terminate at the CP represented by the solid symbols. 

Fig 2 (left panel) displays the spinodal lines, also in the T—fi plane, for the same values of the magnetic field showing 
that very high fields increase significantly the size of the metastable regions. In the T — pB plane, the coexistence 
boundary is limited by the two distinct, low {p^) and high (p^), densities with < p^ for T < Tc and p^ = p^ 
at T = Tc- The boundaries, for different values of B, are shown in Fig 2 (right panel) where we work in terms of 
the baryon density, pg = p/3, given in units of the nuclear matter density, pq = 0.17 fm~'^. An interesting feature 
of this figure is the oscillation of the coexistence boundary around its value at -B = 0, which is more pronounced in 
the p^ branch. The decrease in p^ for eB = 6mJ, at low temperature, shown in the right panel of Fig 2 can be 
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understood in terms of the filling of the Landau levels and, with this aim, we present Fig 3 (left panel) which displays 
the baryonic density and the effective quark mass as functions of the magnetic field at T = 0. To analyze the figure 
let us recall that, in the limit T — ^ 0, the baryonic density can be written^ as [lO] 



PBit^,B)=eikl)J2 Ec^J-^^k^ , (2.12) 



where kp = ^//i^ - 2\qf\kB - A'P and 



or the nearest integer. 



67r2 

f=u k=0 



,„^,,^>l^, (2.13) 

},max 2\qf\B ' ^ ' 




T [MeV] T [MeV] 



FIG. 5: Left panel: Interaction measure for the NJL model as a function of T, for /i = 0, displaying the crossover behavior. 
Right panel: Interaction measure for the NJL model as a function of T, = /ic(B = 0) = 328 MeV displaying the first order 
transition behavior. 
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FIG. 6: Critical regions for the NJL model at eB = (left panel) and eB = 15 (right panel) given by the contour lines of 
the quark number susceptibility, x?, in units of A^. The control parameters have been normalized by Tc and fic corresponding 
to each value of the magnetic field. In both cases the critical region is elongated in the direction parallel to the first order 
transition line. 



^ There is a misprint in Eq. (30) of Ref. |20|| where it should be ps instead of p. 
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Eq. (I2.12[) shows that if < then ps = which is precisely the low density value at T = which is easy to 
understand by recalling that the effective mass is double valued when the first order transition occurs presenting a 
high (M^) and a low (M^) value with < M" for T < and = at T = T^. Now, at T = 0, M" 
corresponds to the value effective quark mass acquires when T — Q and ii — (the vacuum mass) which corresponds to 
~ 403 MeV at B = 0, ~ 416 MeV at = 6 , and M" ~ 467 MeV at eB = 15 . On the other hand, at 
T = the first order transition happens when ^ ~ 383 MeV for B = 0, ~ 370 MeV for eS = 6 and ^ ~ 339 MeV 
for eB = 15m^ so that = even at the lowest Landau level, as required by 9{kp) in Eq. (|2.12p . Then, to 
understand the oscillations let us concentrate on the branch which is shown, together with AI^ (the in-medium 
mass), in Fig. 3 (left panel) where it is clear that both quantities have an opposite oscillatory behavior. The origin of 
the oscillations in these quantities can be traced back to the fact that k^ax (the upper Landau level filled) decreases 
as the magnetic field increases. The first and second peaks, of the curve, correspond to the change from k^ax = 1 
to kmax = for the up and down quark, respectively. As we have seen, for very low temperatures the value of /i 
at coexistence decreases with B (see Fig 1) so that, generally, kmax and M must vary and when k^ax decreases, M 
increases. It then follows, from Eq. p.l2p . that pB must decrease. When kmax = for both quark flavors there are no 
further changes in the upper Landau level and the low temperature oscillations stop at eB ^ 9.5 m^. The magnetic 
field seems to cause a big change on the shape of the coexistence region so that the same high density, pg, may coexist 
with distinct low densities, pg, for two different values of the temperature. For example, p^ ~ 2.5 po coexists with 
pj^ at T Ri 20MeV and with « 1.8 po at T « 70MeV. This pattern is also observed, although in a milder way, 
at the high value eB — 15 to^. Our results also suggest that the highest density is achieved at temperatures close to 
Tc, in opposition to the B = case where this happens at T = 0. Another interesting quantity to be investigated in 
connection with the density oscillation is the latent heat, A£, which at T = is simply given by A£ = p{p^ — p^) 
since the the two coexisting densities occur at the same pressure and chemical potential. Figure 3 (right panel) 
shows that the latent heat value also oscillates around its i? = value for eB ^ 9.5 m^. Figure 4 presents the more 
intuitive P — T phase diagram showing the region of broken symmetry, which corresponds to the "gas" phase, and the 
region of (partially) restored chiral symmetry, which corresponds to the "liquid" phase in analogy with a liquid-gas 
transition. The figure shows that the transition from the "gas" (CSB) phase to the "liquid" (GSR) phase occurs at 
lower pressures when B increases, as expected. Let us now look at some other thermodynamical quantities, such as 
interaction measure (or trace anomaly), which is defined by 

The effect of B over this quantity is shown in Fig. 5 for /i = (crossover) and the critical chemical potential value 
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FIG. 7: The entropy per baryon number, S/A, for the NJL model at _B = (left panel) and eB = 15 (right panel). 
The thin continuous curves correspond to S/A = 5, 10, 15, 20, 40, 100. The thick continuous lines represent first order phase 
transitions while the dashed lines represent crossovers. The focusing effect is not observed. 

at B ~ Q, p,c = 328 MeV (first order transition). The presence of a magnetic field enhances sharper transitions. Next, 
in order to verify how the magnetic field affects the critical region of the phase diagram, let us consider the quark 
number susceptibility, which is defined as 

dp 
dp 

Figure 6 shows this quantity, in units of A^, for S = and eB = 15 m^. In both cases one observes that the critical 
region set by the contour of Xq has an elongated shape due to the enhancement of this quantity along the first order 



Xq 



(2.15) 
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transition line. As discussed by Schaefer and Wambach [43|, at i? = 0, the reason for this behavior is that along the 
path asymptotically parallel to the first order line the quark number susceptibility scales with the exponent 7^ which, 
within the MFA, is expected to be 7^ = f while for any other path the divergence scales with e, which in the MFA 
takes the value e = 2/3 < 7^. Although we do not perform the explicit evaluation of the associated critical exponents 
here we can, nevertheless, take this explanation as being valid also at B ^ since the shape of the critical region 
remains elongated along the first order transition line when a strong magnetic field is present as shown by our results. 
Finally, another interesting physical quantity, that is easily obtained from the EoS, is the entropy per baryon number 

At B ~ 0, the quantity S/A was considered in Ref. [Sjj in order to check an eventual convergence of the adiabats 



toward the critical point as claimed in Ref. [4j|. However, the authors of Ref. [3J| did not observe such an effect 
and the result was explained by recalling that since there is no change in the degrees of freedom of the two phases 
one should not expect the focusing effect to arise within the LSM and the NJL type of models. In this case the 
adiabats show the typical behavior of an ordinary liquid-gas phase transition represented, respectively, by the chirally 
symmetric and broken phases [4^. Our results, displayed in Fig. 7, show that the magnetic field does not produce 
any noticeable effect on the behavior of the adiabats which present a similar pattern at vanishing and high B. 



III. THERMODYNAMICS OF THE LSM UNDER A MAGNETIC FIELD 

In this section, we follow Ref. [2^] and consider the LSM in the approximation where the mesonic sector contributes 
only at the classical (tree) level while one loop quantum corrections, including vacuum contributions, are furnished 
by the quark sector alone. We will perform our analysis considering two sets of parameter values, which produce a 
high and a low value for mo-, in order to address the efficiency of the MFA in dealing with this particular model in 
the presence of magnetic fields. 



A. The Model 



The lagrangian density of the LSM with quarks reads 

>Clsm = V4«^ - .9(^ + «75T • 7?)]^ + ^{d^^ad^a + 9^7? • ^'^tt) - U{a, ir) , (3.1) 

where "0 is the flavor isodoublet spinor representing the quarks, and 

[/(a, Tf) = ^((7^ +Tr2-w2) 2 (3.2) 

is the classical potential energy density. In the chiral limit (obtained by setting /i = in the previous equation) 
the chiral symmetry SU{2)v x SU{2)a is spontaneously broken at the classical level, and the pion is the associated 
massless Goldston boson. Here, we are interested in the h case which implies that chiral symmetry is explicitely 
broken giving the pion a finite mass at T = and /i = 0. The parameters are usually chosen so that chiral symmetry 
is spontaneously broken in the vacuum and the expectation values of the meson fields are (tr) = /,r and (tt) — 
where f-^ — 93MeV is the pion decay constant. When pure vacuum contributions are considered within the MS 
renormalization scheme, as we do here, an arbitrary mass scale (A^) also appears in the final results. As discussed 
in the introduction, for our purposes, it will be necessary to consider two different sets of parameters in order to get 
a high and a low value for TOq-- The first set is given by Aj^ — 16.48 MeV which, together with v = 64.29 MeV, 
A = 46.06, and h = 1.77 x 10^ MeV^ , yields = 138 MeV and = 600 MeV. The second set is just like the first 
except for A and v which are set to A = 36.96 and v — 54.96 MeV yielding 111^ = 450 MeV. In both cases the sigma 
meson mass falls within the range of the broad resonance nia = 400 — 800 MeV [1^ . 



B. The LSM free energy in the presence of a magnetic field 

Treating the bosonic degrees of freedom at the tree level implies that the loop contributions to the free energy come 
entirely from the fermionic sector. This contribution is represented by an integral similar to the one which appears 
in Eq. (|2.2[) so that the free energy is given by 
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FIG. 8: The LSM phase diagram in the T — fi plane for B — 0, eB — 6m^ and eB — 15 m^. The left panel, in which 
only crossover occurs, refers to rua ~ 600 MeV and the right panel to rricr = 450 MeV. In the right panel, the crossover 
is represented by dashed lines, the first order phase transitions by continuous lines and the CP, by solid symbols which 
are located at (T^ ~ 79.32 MeV, ^ 198.52 MeV) for B = 0, {T, ~ 12.3 MeV, 343.4 MeV) for eB = 6m^, and 
(Tc ~ 34.2 MeV, /ic — 365.1 MeV) for eB = 15 m^. In both cases the region of broken symmetry expands with increasing B. 
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FIG. 9: The quantity —d{a)/dT, where {a) is the LSM order parameter, as a function of T for rria = 450 MeV. The left panel 
corresponds to n = and the right panel to fic{B = 0) = 198.52 MeV. 



7-LSM rr/ -;\ , 7- 



LSM 
vac 



LSM 
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(3.3) 



where, the magnetic and the in-medium terms bear the same form as their NJL counterparts with the obvious 
replacement — > = S^ip^ +7?^). On the other hand, the pure vacuum contribution, whose general from is 
given by Eq. (|2.7|) . is treated in different fashion within both models and in many LSM applications it is discarded 
on the grounds that it does not give significant contributions at high temperature and/or chemical potential values. 
However, a considerable number of studies have pointed out to its importance in the characterization of the phase 
transition Here, we follow Ref. [28| by including the pure fermionic vacuum contribution in the free 

energy. Since this model is renormalizable the usual procedure is to regularize divergent integrals using dimensional 
regularization and then to subtract the ultra violet divergences in the MS renormalization scheme. This procedure 
gives the following finite result 



LSM 



(8^)2 



In 



. MS, 



(3.4) 



The expectation values, (cr) and (tt) fields are obtained by extremizing the free energj 



ax-LSM ax-LSM 

, = 



da 



(3.5) 
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FIG. 10: The sigma meson mass, 771^, whose inverse represents the correlation length, as a function of T for = (left panel) 
and ^c{B = 0) = 198.52 MeV (right panel). For comparison reasons the mass has been normalized by its value at T = 
(mo- = 450 MeV). The curve correspondig to B = in the right panel shows that rria — >■ when /i — > /ic(B = 0), as expected. 




FIG. 11: Left panel: Manifestation of magnetic catalysis, at T = 0, for the NJL model (continuous line) and for the LSM 
(dashed line) with mo- = 450 MeV. The effective quark mass, Mq{B), which represents the effective quark mass for each model, 
has been normalized by its B = value to allow for comparisons. Right panel: The pseudocritical temperature, as a function 
of B, for each model. This quantity has also been normalized by its value at B = for comparison reasons and the LSM is 
taken with nia = 450 MeV. 



while the mesonic masses are determined by the curvature of the free energy at the global minimum 



2 _ a2_^LSM ^ a^^LSM 



evaluated at cr = (ct) and tt = (tt). Here, we follow the usual procedure [EfiUll by setting (t?) = so that — g{a). 



C. Numerical Results for the LSM 



Let us start by mapping the T — /i phase diagram, at the physical point, for three relevant values of the magnetic 
field {eB = 0,6m J, 15m^). Let us start with the parameter set which gives m„ = 600 MeV whose result is shown 
in the left panel of Fig. 8 which indicates that the high-T/low-/i behavior is similar to the one found in the NJL 
model. Namely, a crossover happens at Tpc values which increase with the magnetic field. Note that, qualitatively, 
this situation does not change even by going beyond the MFA or employing other parametrizations |29|, Is^]- Here, 
the results given by the parametrization leading to rria = 600 MeV start to depart from the ones observed in the NJL 
case at intermediate values of T and fj, since the crossover does not change into the expected first order transition 
with its associated CP. Then, the whole T — ^ plane is dominated by the crossover exactly as the authors of Ref. 
[28t have observed using rric = 800 MeV. The situation changes drastically if one uses the parameter set which yields 
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nia = 450 MeV as the right panel of Fig. 8 suggests. Indeed, with this set one observes the crossover/CP/first order 
transition pattern. However, when the magnetic field increases one does not observe the characteristic effects observed 
earlier in the NJL model as the decrease of the coexistence chemical potential at low T (IMC), the increase of Tc and 
the decrease of /ic for example. We can further investigate the strength of the transition by analyzing how the order 
parameter varies with the temperature as shown in Fig. 9 for = and fJ,c{B = 0) = 198.52 MeV. In both cases one 
sees that the strength, characterized by the peak of —d{a) /dT decreases with increasing magnetic fields and the figure 
also suggests that, when B ^ and n ^ 0, the crossover becomes smoother as B increases. Based on the findings of 
Ref. [23], which show that a magnetic field makes the transition sharper, one would expect the opposite behavior. 
Since the correlation length is governed by 1/mo- [lH the sigma meson mass, as a function oi B,T and fi is also an 
interesting quantity to be analyzed and with this aim we present Fig. 10. In Ref. [sOl, which considers the PQM at 
/i = 0, it has been demonstrated that the value of rrirj at Tpc is almost independent of B (both, in the MFA and FRG 
approaches). Contrary to that, within our approximation, we observe that the difference ^^(Tpc) — TOo-(O) decreases 
with increasing B which is also in opposition to the expectations furnished by the results of Ref. [l^l- Within the 
NJL model the behavior of m,cr{T, fi, B) has also been recently addressed [1^1 furnishing results which predict that 
mcr(Tpc) — nicr(O) increases with increasing B and also that the transition appears to be sharper as /i increases. These 
comparisons suggest that the naive application of the MFA to the LSM at finite T, fi and B indeed misses some 
important physics, especially at finite densities. This conclusion is not only supported by the NJL-MFA results of 
Refs. [22 - l2a | but also by the LSM-FRG results obtained by Andersen and Tranberg |29i|. 



D. Comparing Magnetic Catalysis in the LSM and the NJL model at = 

In the previous subsection we have pointed out that the MFA does not seem to furnish reliable results for the LSM. 
However, this issue seems to be less severe at = where at least the predicted type of transition is in agreement 
with other model applications as well as with lattice simulations. We can then address the question of how magnetic 
catalysis [l^] (see Ref. [1^ for an updated review) affects each model by comparing the dimensionless quantities 
Tj^/Tpj. and Mq{B)/Mq{0) where, now, Mq is representing the effective quark mass in each case, that is, M in the 
NJL model and in the LSM. The result, shown in Fig. [TTl suggests that within our approximations the effective 
quark mass of the LSM is more sensitive to the presence of a magnetic field. At eB — 15 the LSM quark mass 
has increased by about 20% while the NJL quark mass has increased by about 15%. The increase of Tpc is even more 
dramatic within the LSM and at eB — 15 it increases by about 50% in relation to the B = value while in the 
NJL the increase is about 9%. At low values of B the LSM pseudocritical temperature suffers a sudden increase which 
again may be an indication of the failure of the approximations we have adopted for this particular model. 



IV. CONCLUSIONS 



We have considered both, the LSM and the NJL model, in the description of hot and dense two flavor quark 
matter subject to strong magnetic fields, such as the ones expected to be created in non central heavy ion collisions. 
Here, the simplest version of these models, which do not consider the Polyakov loop, have been treated within the 
MFA. In this approximation, one loop fermionic contributions have been considered, allowing for the determination 
of the phase diagram associated with chiral symmetry breaking/restoration. So far, most applications [H, [H, [H, [2^ 
at finite T, /i and B have mainly focused in the determination of the T — ii phase diagram without analyzing the 
effects of the magnetic fields over important physical quantities such as the coexistence boundaries, the quark number 
susceptibilities, the interaction measure and the mesonic masses among others. Therefore, our aim was to perform 
this analysis in order to understand the effects of B in more detail, especially at low temperatures, which is the less 
explored region. We have started by considering the NJL model, with physical quark masses, obtaining a T — phase 
diagram which is in qualitative agreement with the one obtained in Ref. [2^ for three quark flavors and also with 
Ref. [26| for two quark flavors. One of the most important features of this diagram is related to the position of the 
critical end point that appears to be located at higher temperatures and smaller chemical potentials when B increases. 
One also observes that the first order segment of the transition line increases with B while the coexistence chemical 
potential value decreases at low temperatures. These observations suggest that the magnetic field has a direct effect in 
the physics associated with the first order transition and may influence the size and shape of the coexistence region for 
example. In order to check that we have mapped the T — fi plane into the T — pB plane showing that for eB < 9.5 
the high density branch of the coexistence phase diagram oscillates around its i? = value as a consequence of filling 
the Landau levels which influences the values of quantities such as the latent heat as we have shown. This finding 
may also have consequences regarding, e.g., the physics of phase conversion whose dynamics requires the knowledge 
of the EoS inside the coexistence region. For example, at a given temperature, the surface tension between the two 
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coexisting bulk phases {p^ and p^) depends on the value of their difference [i^, [13] which will, possibly, be affected 
by the oscillations suffered by the coexistence boundary due to the presence of a magnetic field. This feature of 
the phase diagram deserves further investigation since it may be important in the description of the EoS inside the 
boundaries. Then, using the quark number susceptibility, we have mapped the critical region observing that its shape, 
which is elongated along the first order transition line, is not affected by the presence of a magnetic field. We have 
also investigated the behavior of the adiabats in the presence of B observing the same pattern which was observed 
at i? = [3J1. The two latter results indicate that general model characteristics, associated with the shape of the 
critical region and the isentropic trajectories, are unaffected by the presence of B. For the NJL model we conclude 
that magnetic fields appear to cause a significant effect at the intermediate and low temperature parts of the phase 
diagram causing the crossing of the transition lines representing different values of B. The presence of magnetic fields 
also promotes a change on the size and location of the first order line. 

Our results for the LSM indicate that the MFA treatment as performed here does not seem to be appropriate, 
especially at higher values of the chemical potential where, as discussed, the magnetic field has a strong influence. 
In the light of the results furnished by the more powerful FRG [H, 1131 and also those obtained with the NJL model 
within the MFA [l^-HBl we have observed a significant disagreement, which include the absence of IMC in the LSM 
phase diagram, in the low temperature region. 
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